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Abstract. 

Recent results indicate that the J /if/ suppression pattern differs with rapidity showing a larger 
suppression at forward rapidity. J /if/ suppression mechanisms based on energy density (such as color 
screening, interaction with co-movers, etc.) predict the opposite trend. On the other hand, it is expected 
that more cc pairs should be available to form quarkonia at mid-rapidity via recombination. Some models 
provide a way to differentiate J/if/ production from initially produced cc pairs and final state recombination 
of uncorrelated pairs, via the rapidity and transverse momentum dependence of the elliptic flow (V2). 



During 2007 data taking at RHIC, a large sample of Au+Au collisions at -\/s NN =200 GeV was collected. 
The statistics has been increased compared to previous 2004 data set, thus allowing a more precise 
measurement of the J /if/ production at both mid and forward rapidity. Furthermore, the PHENIX 
experiment benefited from the addition of a new detector, which improves the reaction plane resolution 
and allows us to measure the J /if/ v%. Comparing this measurement to the positive D-mesons v 2 (through 
non-photonic electron decays) will help constraining the J/ if/ production mechanisms and getting a more 
precise picture of the proportion of J/if/ coming from direct production or charm quark coalescence. 
Details on how the J/if/ V2 is measured at mid-rapidity rapidities are presented. The J/if/ V2 as a function of 
transverse momentum are compared to existing models. 



1. Introduction 

The use of charm quarkonia produced in ultra-relativistic heavy-ion collisions as quark gluon plasma 
(QGP) probes has a rich tradition, starting with the seminal paper by Matsui and Satz HI, which predicted 
the suppression of J/i// production in a QGP. At RHIC, PHENIX has measured the J/if/ production in 
A-A collisions at t/saw = 200 GeV. The suppression as a function of the number of participants is 
found to be similar to the one at lower energy at SPS [2] even though the energy in the center of mass 
differs by an order of magnitude, leading to different cold nuclear effects (CNM) as well as the melting 
of other excited charmonium states. On the other hand, the larger energy and charm quark densities 
created at RHIC may allow regeneration of the J /if/ via recombination of uncorrelated charm quarks 
during the collision 0IHE]. Moreover, PHENIX recently published results [6] indicate that the nuclear 
modification factor Raa is smaller at forward rapidity than at mid-rapidity. This may come from the 
amount of CNM effects that may be different for each rapidity [7 ], or from coalescence of charm quarks 
which would enhance the J ft// production further at mid-rapidity. 

Studying other observables, such as the azimuthal anisotropy of the produced particles might give 
additional information on how these particles are produced [8 ]. In non-central heavy ion collisions, the 
distribution of the colliding matter is not isotropic around the beam direction. In the thermodynamic 



picture, the asymmetric distribution of initial energy density causes a pressure gradient which is larger 
in the shortest direction of the ellipsoidal medium. Therefore, the momentum distribution of produced 
particles can also be anisotropic. PHENIX measured a positive elliptic flow for heavy (c,b) quarks via 
the study of so called non photonic electrons for a transverse momentum higher than 1 GeV/c 0. If 
J /tf/ are produced from cc recombination in a deconfined phase [ 10], they should inherit their flow. On 
the contrary, J/</r directly produced by hard QCD processes early in the collision cannot be sensitive 
to collective phenomena. Therefore measuring a positive J /iff elliptic flow would indicate the level of 
recombination that takes part in the J /iff production mechanisms iTTTTl . 



2. J /iff elliptic flow 

The elliptic flow is quantified by the second Fourier coefficient v 2 of the J /iff azimuthal angle 
distribution lfT2l[T3l : 

* N „ = M + 2 • v 2 • cos(2 ■ (cf> - iff))] (1) 

with <p the azimuthal angle of the particle with respect to an arbitrary axis (Oj), and if/, the reaction plane 
angle (plane defined by the impact parameter vector and the beam axis) with respect to (Oj). 

The reaction plane of each event is measured using the azimuthal distribution of a limited subset 
of particles produced in the collision. Its resolution depends on the number of particles detected, the 
strength of their flow, and the resolution of the detector. Therefore, the elliptic flow measurement is 
corrected so that this resolution is reflected on the true v 2 measurement: v'^ ue = v™ eas /ctrp. 

Starting from 2007, PHENIX uses two new Reaction Plane Detectors (RxnP) (see Fig.Q]) to measure 
the reaction plane of each collision following methods described in lfT"3l [T4l . It improves the reaction 
plane resolution, and thus gives a correction ctrp = (cos(2 Aif/Rp)) twice better than what was achieved 
previously with the Beam Beam Counters (BBC), or what can be measured with the Muon Piston 
Calorimeters (MPC) (see Fig. |2]> . 
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Figure 1. PHENIX detector during 2007 data taking 
with the RxnP detector near the collision vertex. 
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Figure 2. Reaction plane resolution as a 
function of centrality, measured with the RxnP 
(top), the MPC (middle), or the BBC (squares). 



The PHENIX experiment is designed to detect heavy quarkonia at forward rapidity (1.2 < \y\ < 2.2) 
via the yu + yu~ decay channel and at mid-rapidity (\y\ < 0.35) via the e + e~ decay channel ifTBl . At forward 
rapidity, cathode strip tracking chambers and alternating layers of steel absorber and Iarocci tubes allow 
muons to be identified and tracked over an acceptance of A<p = 360°. At mid-rapidity drift chambers, 
ring imaging Cerenkov detectors and electromagnetic calorimeters are used to detect electrons in two 



arms covering A<p = 90° each. During the 2007 run at the Relativistic Heavy Ion Collider, PHENIX 
recorded Au+Au collisions and 1.8 x 10 9 events are analysed here (42% of the total luminosity). 

The forward measurement of J /if/ elliptic flow benefits from three times more J /if/ than at mid-rapidity 
but with a poorer resolution (150 MeV vs. 80 MeV) and signal over background ratio (1/6 vs. 1/3). In 
addition, at mid-rapidity both RxnP detectors can be used leading to an optimal measurement of v%. A 
challenge to the forward rapidity measurement comes from the overlapping acceptances of the RxnP with 
the muon arms (see Fig.Q]). The two muons and all possible accompanying fragmentation products might 
go through the RxnP, which may introduce non-flow effects and auto-correlations between the reaction 
plane and the J /if/ azimuthal angle. To limit this bias on the V2 measurement, the best configuration 
uses the RxnP located in the muon arm opposite to the J /if/. The drawback is that the resolution of the 
reaction plane measurement is worse than when using both forward and backward RxnP detectors. On 
the other hand, using both RxnP or restricting to the one where the muons go (and thus possibly biasing 
the measurement) could be used to estimate the magnitude of the bias. 

J/if/V2 results are still underway at forward rapidity even though the Raa has been measured lPT6ll and 
is compatible with published results Q. The method consists in extracting the number of J /if/ counts 
in bins of <p-if/, transverse momentum pj and centrality using mixed event background subtraction. The 
v meas ^ s extracted from the fit of the invariant yield with Eq. Q] in each pj and centrality bin, and is 
corrected with cr RP . A similar method was also applied for mid-rapidity, but the result presented here 
uses an alternate method based on the additivity of the elliptic flow ffTTl , which provides a more precise 
handling of the background elliptic flow: 

v S 7 (Mjm) - < G (M JM ) • „ W - / BG (M JM ) W (2) 

v fg (Mj/,i,) is the V2 of the foreground, from same event unlike sign pairs, in [2.9, 3.3] GeV/c 2 . A^ FG (Mj/^,) 
is the number of counts in the 3 /if/ mass range from mixed event pairs. N s {M^i^) is the number 
of measured J/if/, N s = N FG - (N comb + N rem ), with N comh the contribution from the combinatorial 
background estimated using event mixing, and N rem the remaining contribution of the background 
(physical processes or uncertainty on the mixed background normalization) fitted with an exponential. 

f BG (M) is fitted outside the J /if/ mass to v FG (M) ■ ^S^ , and extrapolated to [2.9, 3.3] GeV/c 2 . 



3. Results 

The J /if/ \>2 at mid-rapidity for 42% of the Run-7 statistics is shown for three pj bins in the centrality 
range [20,60%] on Fig. [3] Statistical errors and point to point uncorrected systematic errors coming 
from the signal extraction are drawn with vertical bars. Systematic errors correlated between the pj bins 
appear with boxes. These errors account for the V2 extraction and are estimated by the difference on the 
V2 value when fitting with Eq.[T]or using directly {cos2(<p - iff))- They also account for the background V2 
shape, estimated with a fit to the background distributions with a constant, linear or quadratic polynomial. 
Global systematic errors to account for the technique used to determine the reaction plane angle and 
resolution are written on the figure. 

The measurements shown on Fig. [3] have big uncertainties mostly due to lack of statistics and 
to the poor signal/background ratio. The V2 integrated over all pj is compatible with 0: V2 = 
-0.10 + 0.10 ± 0.02(±3%). Some theoretical scenarios appear on Fig.[3]to illustrate different production 
mechanisms predictions. When only initially produced J /if/ (no regeneration) are considered lTT8l with 
normal suppression from inelastic collisions with spectators and anomalous suppression from collisions 
with gluons in a QGP, the prediction is slightly positive but reaches at most 0.02 for pj = 5 GeV/c. 
Models with partial or full charm quark or charmonium thermalization should have larger V2 than this. 
The addition of recombination of c and c during or at the term of the partonic phase should enhance 
the charmonium production and affect the elliptic flow. For instance, full regeneration of the J /if/ at 
freeze out is studied in [19] and in a transport model description in ll20l . Those models predict a steep 
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Figure 3. J/i// v 2 at \y\ < 0.35 for [20,60%] in centrality, as a function of p T , using 42% of 2007 Au+Au 
statistics, with some predictions. Error bars are statistical errors uncorrected point to point, boxes are 
systematic errors correlated between each pr bins, and a global systematic error is written (see text). 



increase of the J/t// v 2 , reaching 0.2 at pj = 5 GeV/c. If the production is a mixture of regeneration and 
initial production as in lF2Tl l22l . the predicted v 2 has a more moderate slope and reaches intermediate 
values. Finally, comover interactions |[23l leave the J/i^s with no vi within errors, slightly positive at 
pj = 2.5 GeV/c, and negative at pj = 4 GeV/c. The mid-rapidity J/i/f v 2 measurement obtained with 
less than half of the potential 2007 statistics does not allow to distinguish between predictions. Work is 
ongoing to process the remaining statistics and to reduce the systematic uncertainties at mid and forward 
rapidities. Future RHIC Au+Au runs will provide more statistics which is crucial for this measurement. 
The experimental study of J/i/f v 2 will be complementary to other studies of J/i/f suppression and will 
provide valuable information about the early stage of high-energy heavy-ion collisions. 
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